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Abstract 

We study the free energy of the 1 + 1 dimensional O(N) nonlinear a- 
models for even N using the TBA equations proposed recently. We give 
explicit formulae for the constant solution of the TBA equations (Y-system) 
and calculate the first two virial coefficients. The free energy is also compared 
to the leading large iV result. 



1 Introduction 



A large class of integrable two-dimensional models can be described as perturbed 
conformal field theories (PCFT). If the perturbing relevant operator is chosen so that 
perturbation preserves integrability then the spectrum of the resulting (in general 
massive) integrable model and its exact S-matrix can be conjectured Q. Further, 
the set of TBA equations describing the ground state energy in finite volume (or 
equivalently the free energy of the system at finite temperature T) can be written 
down. 

Non-linear cr-models (NLS) are another important class of integrable models. 
Their spectrum and S-matrix have been obtained by bootstrap methods |5| but the 
corresponding TBA system was not known until recently. In ||, |[ (based on earlier 
constructions [7], ^|) NLS models are represented as limits of PCFT models. More 
concretely, the NLS model whose fields take value in the coset space G/H can be 
described as the k — > oo limit of the (G/H) k conformal coset model || perturbed by 
a certain relevant operator. This way of representing NLS models gives additional 
insight into the structure of the model and is reminescent of the description of the 
model as the infinite flavour limit of a certain fermionic model ||1U|| . Moreover, it 
allows one to calculate the ground state energy of the NLS model by considering 
the limiting case of the TBA system, available for finite values of the Kac-Moody 
level k. 

Because of the obvious importance of the NLS models in this paper we perform a 
number of tests for the above TBA system. For simplicity, we restrict our attention 
here to the O(N) NLS models for even N. We show that in the T — > oo limit 
the correct UV central charge is reproduced. In the opposite (low temperature 
or equivalently low density) limit there exists a systematic expansion, the virial 
expansion. It is known |jTT| that the virial coefficients are completely determined 
by the scattering data alone. This is useful since the virial coefficients can also 
be calculated from the TBA equations by solving the linearized system in Fourier 
space. We work out the first two virial coefficients directly from the S-matrix and 
also from the TBA system and find that they agree completely. Finally, NLS models 
can also be solved in the 1/N expansion. We calculate the free energy in the leading 
large N limit and compare it to the N — > oo limit of the virial coefficients calculated 
previously. 

2 TBA equations 

Here we study the TBA system proposed in |4[] for the O(N) nonlinear sigma model 
for the case of even N. The TBA equations are formulated in terms of the densities 
Eq(6) and £$(9), where the index a labels the nodes of the D r Dynkin diagram 
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(N/2 = r > 2) and the range of the lower index is m = 1, 2, . . . , k — 1 (k > 2). Let 
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where (fl|) and (^) are defined for all a = 1, 2, ... r and m = 1, 2, . . . , k — 1. By 
convention we take B { a) (6) = B { k a) (6) = 0. g = 2r - 2 is the Coxeter number of D r . 
We further define 
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where for r > 3 AC^fu;) is the matrix inverse of 
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Here I a b is the incidence matrix of the D r Dynkin diagram. 
For later use we note that N al (u) is given by 

a , 2r-2-a 

Ki(u) = \ \„_ o = l,2,...,r-2 
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and 
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where q = e ~ . We will assume that (0) holds also for r = 2, by definition in this 
case. 

The TBA equations proposed in H are 
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and for a = 1, 2, . . . , r and m — 1, 2, . . . k — 1 
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In (||) M is the mass of the particles and if is a chemical potential coupled to particle 
number. 

Finally the pressure of this one-dimensional gas is given by 

p = / d9 cosh In [1 + e- £oW ] . (10) 

The above TBA system describes the perturbed (0(N)/0(N — 1))& coset model 
for finite k. The pressure of the O(N) NLS model can be obtained from it by taking 
the k — > oo limit Q . 



3 Virial expansion 

We have introduced the chemical potential H to facilitate the virial expansion of 
the pressure. In the limit H — > — oo the densities can be expanded in powers of the 
fugacity y = e u l T as 

1 oo 

eo(0) = -(Mcosh£-tf) + ^4^(% M , (11) 
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whereas the pressure is expanded as 

oo 

p = tJ2p (>1 V- (13) 
11=1 

Alternatively, one can get the same expansion coefficients (without introducing the 
chemical potential) by considering the low temperature expansion of the pressure. 
In this case the coefficient is of order exp(— /iM/T). 

It is known |11[ that for the case of a system of iV particles (independently of 
integrability) the virial expansion coefficients are determined by the S-matrix. If all 
particles have the same mass M then the first two virial coefficients are 



j9 (1) = — / d#Mcosh#e-^ coshe (14) 
and p^ = p± + p^ , where 



p ( 2 ) = _^/ d£Mcoshfle-^ coshe (15) 
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and 
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Here n = N, <pi(9) = -i^lnSj(0), where {^(0); / = 1, 2, ... , iV 2 } are the eigen- 
values of the 2-particle S-matrix and s is a factor related to particle statistics. For 
the case of diagonal scattering, s = =Fl for bosonic/fermionic type particles |jT2| , 
whereas for all known non-diagonal cases (where all particles are of fermionic type) 
s = N. Finally, for later use, we introduce the Fourier transform 

/oo 
due+*}(u>). (17) 
-oo 

For the O(N) nonlinear sigma model, using the eigenvalues of the bootstrap S-matrix 
we get 

<f>(u) = N e" wM - —— - e"^. (18) 
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4 Solution of the Y-system 

Let us consider the generic r > 3 case first. The leading terms in the expansion ([12] 
turn out to be ^-independent and the constants 



i(o) 



can be determined by solving the equations 
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for a = 1, 2, . . . , r, m = 1, 2, 



fc — 1, where by convention x Q 
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ubiquitous equation is called the constant Y-system and its solution is well known for 
all classical Dynkin diagrams |L3|. It is usually solved by relating it to the Q-system 
equations 



m+l Qm-l 



6=1 



for a — 1, 2, . . . , r, m — 1, 2, . . . , k — 1, with "boundary condition" 

Q< a) = Q[ a) = l. 
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We are interested in positive solutions of fl2~T|), related to (|20| ) by 
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We will need the explicit solution of fl2"ip only for a = 1. This can be written as 
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(); u = , ;;;ffi / 2r " ;/ and o 



(24) 



(i) _ p(r)p(2r - 2) ^ (1) _ p(r + l)p(2r - l)p(2r - 2) 

p(l)p(r — 1) *~ p(l)p(2)p(r — 1) 



(25) 



We can now calculate the first virial coefficient. It is of the form (|1J), where the 
coefficient n can be expressed in terms of the constants (|T9|) as 
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If we use (0), (0) together with 



can be written as 
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From fl25|) we see that lim n = 2r = N. 
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5 The second virial coefficient 

Encouraged by the agreement of the first virial coefficient with the expected form 
we now go one step further and calculate the second virial coefficient. Using the 
expansions ([□]) and ( |I2"D in (|4"5|), ( [46]) and going to Fourier space we find that p^ 
is of the form (|T5| ) with s = n and is of the form ([H]) with 



(w) = n 2 £j\U")S$ B) (?)- 
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Here the functions S[ a \q) can be calculated by solving the linear equations 

(? + Q' 1 ) (l + 7^) ) ^ = 5 ^ml + + + £ I^S® (29) 



for a = 1, 2, . . . , r, m = 1, 2, . . . , k — 1, where 5 
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The solution of the linear equations ( p9|) can be given in terms of the numbers 
Qm satisfying (pT|). Although their explicit expression is known [Tj|, it is actu- 
ally easier to verify the validity of the solution given below by using the Q-system 
equations ( pTD directly. 

The solution for S$ (q) is given as follows. Using the definitions 
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for m = 1, 2, . . . , k — 1 (for m = 1 the second term in (|33|) is absent) together with 



r(a) 



r(a) 



we first define 
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for m = 0, 1, . . . , k. Next we define 
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Note that both Hm (<?) and Sm^ (q) are finite polynomials in q. 
The solution of (ED) is finally given by 
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Using the solution (|38| ) we can calculate 

1 1 



where 



X _|_ q9 l _ g2k+2g 
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Putting everything together we have 

k{q) = Q« g 2 + (n 2 - - (n 2 - Q«) g»+ a - g 2 * 

and using also 



lim = (r + l)(2r - 1) 
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we see that fl39|) indeed reproduces (ITS) in the limit k — ► oo. 



6 Calculation of the UV central charge 

For this purpose it turns out to be useful to rewrite the TBA equations a little. Let 
us introduce the new variables 



eo(0) = e (6), 



L {e)=ln [l + e- eoW ] 
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In terms of these new variables the TBA equations take the standard form 

1 r 
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Our convention for convolution is 
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(f*9) (e) = ^J°°de'f(e-e')g(e') 
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and (6) = L^' (6) = by convention. 
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We note that the physical particles correspond to eo- They are of mass M and 
are coupled to the chemical potential H. All other excitations corresponding to the 
densities are "magnons" . They are massless and are not coupled to the chemical 
potential. 

The standard form of the TBA equations allows us to calculate the UV central 
charge the standard way [ 12 ] . It can be written as cuv — — & 2 ' with 
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where C(z) is Roger's dilogarithm. In ( |4T)| ) the set of numbers {i™} satisfy ( p0| ) and 
are the same as occuring in the virial expansion. They can be found in |13| where 
not only their values, but also the corresponding sum of dilogarithms is evaluated 
(conjectured): 



.(2) 



k(k — l)r 



(50) 
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The numbers {if } are defined as the solution of a similar problem. They satisfy 
the equations 
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for a = 1,2, 
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Similarly to the previous case, it is also useful to introduce the "Q-system" here. 
We denote the corresponding numbers by {Rm^} and define for m = 0, 1, . . . , k 
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The numbers Rj$ are assumed to be nonnegative and satisfy the Q-system equations 

(m^R^R^+mm 1 ^ (55) 

6=1 

for a — 1,2, ... ,r, m — 0,1, ... ,k with boundary condition 



6=1 



We have found the solution of the modified Q-system (p5|), (|56|) . The components 
corresponding to the fork part of the Dynkin diagram can be written as 

• H(fc-mK+p)7r r -2 ■ (k-rn+2l+l)ir . (2/+1)tt 

D(r-l) _ c)(r) TT _ fc+g-1 TT 2(fc+g-l) bm fc+g-1 

-"-m -n-m ~ 1 1 . (a|p)7r 1 1 • (2/+1)tt . (fc-m+2Z+l)7r " \° ' > 

«eA + S111 k+g-1 1=0 S111 2(fc+g-l) S111 fc+g-1 

Here A + is the set of positive roots, p is half the sum of the positive roots, uj r is 
the r th fundamental weight and scalar product is normalized in weight space so that 
(a\a) = 2 for the roots a. For m = the ratio 

sin (fc-^+g- 1 ^ 

. ;- g i): (58) 

SUl t^- 

fc+g-1 



has to be omitted from the first factor in ( p7\ j together with the simultaneous omis- 
sion of its inverse from the last factor. 

The remaining components, {i?™ for a = r — 2, . . . , 2, 1}, can be calculated from 
the equations ([55]) with a = r — 1, . . . , 3, 2. We have verified that the last Q-system 
equation (for a = 1) is then also satisfied. 

Based on extensive numerical study, we conjecture that if we use the solution of 
the "modified Y-system" ([H]), ( |5*2"|) obtained this way in the dilogarithm sum (|S) 
we get 

(1) = fc(fc-l)r + fc(r-l) 

k + 2r-3 K J 

Putting (H) and (H) together, we see that cuv can also be written in the 
suggestive form 

cuv = c(k,N)-c(k,N-l), (60) 

where 



k 



N(N-l) 



is the central charge of the O(N) WZNW conformal model at level k. 
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7 The 0(4) model 



The r = 2 case is somewhat special and we discuss it separately. From (|j) and (|7|) 
it follows that in this case 

We now define the new variables 
e (6) = e„(0), e m (0) = e_ m (^) = -e^(0). (63) 

We further define 

/ m (^) = ln[l + e-^^] (64) 

for m = 1 — k, . . . , 0, . . . ,k — 1 and l±k(0) = 0. This allows us to write the TBA 
equations and @ in the compact form 

e m + i> * (lm+i + lm-i) = ^(Mcosh6» - H)8 m0 (65) 

for m = 1 — k, . . . , 0, . . . , k — 1. 

In this form it is easy to recognize that fl6"5]) is the standard TBA system [0 asso- 
ciated to the A 2 k-i Dynkin diagram with m = corresponding to the single massive 
node. This observation allows us to determine the UV central charge immediately 
from 

c uv = c{A 2k -i) - 2 c(A fc _i). (66) 
Using the well-known formula |L3] 

we get 

3A; 2 

c »v = (fc+W+2) ' (68) 

which is the expected result since it can also be written as 

Cuv = c (k,4) -c(2fc,3). (69) 

Note that the case of the 0(4)/ 0(3) sigma model is exceptional in that the KM 
level corresponding to 0(3) is doubled. For N = 3 ( |61~D is not valid either, we have 
c(k, 3) = 3k/(k + 2) instead. 
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The Y-system equations in this case simplify to 

x 2 m = (1 + x m+ i)(l + x m -i) (70) 
for m = 1,2, . . . ,k — 1 and Xq = x^ = 0. They are solved in terms of the Q-sytem as 
x m = Qm+iQm-i, where Q 2 m = l + Q m+ iQ m _i. (71) 
The Q-system has a simple solution in this case. It is 

p(m + l) 

^ m p(l) ' U j 

It is easy to calculate the virial coefficients and They are of the form (|14"D 
and fllBD respectively with 

s = n =(l+x 1 ) = Q?, (73) 

which indeed reproduces iV = 4 in the k — > oo limit. 

The Fourier transform (p~7|) entering the second virial coefficient in this case is 

2n 2 q 
1 + q 2 

where Si(q) is determined by solving 



^u) = - rT \S l {q), (74) 



q+ - ) ( 1 + — ) 'S'm(g) = 5 m i + 5* m+ i(g) + S^-i^) (75) 



for m = 1, 2, . . . , k — 1 with boundary condition S'o(g) = Sfc(<?) = 0. The solution of 
(HD is 

1 



Sm(<7) = 77- 



Qi l - g 2fc+4 

for m = 1,2, . . . ,k — 1, where 



(76) 



Qm+l m Qm-l m+2 



S m (q) = ^q m -^q m+I . (77) 



Finally we find 



^) = 1 2 "^i-^ [ ft -" 2 + " 2i(1 - fe " 2) ]- (7f 



which reduces to the expected result 



~ 8q 2 (3-q 2 ) 

^ = 1 + q 2 ( ?9 ) 



in the k — > 00 limit. 
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8 Large N limit 

Finally, we consider the N — > oo limit, where the O(N) a-models are exactly solv- 
able. This limit was studied in |14|], where the exact formula 

In ^ = F(z) = 2 £ K (mz) + O (1) (80) 

m=l ^ ' 

for the mass gap M(R) in finite volume (corresponding to a circle of radius i?) 
was found. Here K Q is a modified Bessel function and z = M(R)R, which can be 
obtained, as a function of ( = MR, by solving this implicit equation. 

Using similar techniques, we calculated the N — > oo limit of the pressure. This 



requires the calculation of the ground state energy in finite volume [0. We find 



AT AT 00 

p(R) = — {M 2 — M\R)} + -M\R) V K 2 (mz) + 0(1), (81) 

m=l 

where -K2 is a modified Bessel function and the temperature is given by T = 1/R. 
We first study the R — > limit of the effective central charge 

c( J R) = —p(R). (82) 

7T 



For short distances F(z) ~ 7r/z implying [|1J] 



~ — + (83) 
In (I 

consistently with AF perturbation theory. This can be used to determine the short 
distance (or high temperature) limit of the effective central charge. We find 




c{R) = N{l-—j-^ + 0\ — T ^ | }+0(l), (8J) 



as expected. 

In the opposite (R — > oo) limit flSOD can be expanded as 



= i + X) MO. ( 85 ) 



M 



where A„(£) = (9(e The first two terms of this expansion are 

A X (C) = 2K (C) and A 2 (C) = 2K 2 (() + 2K (2() - 4Ctf (C)#i(0- (86) 
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Similarly, the pressure fl8~T|) can be expanded as 



P(«) = ^EP M (0, (87) 
^=1 



where p^(C) — 0(e vC -). As explained above, these expansion coefficients must be 



the same as the virial coefficients of (|T3|). Using fl36|) in the expansion of (81), for 
the first two expansion coefficients we get 

(i) (c) = ^^ l(C) + 0(1); (88) 

7T 

P (2) (C) = ^^ 1 (2C)-2Cir 2 (C)]+O(l). (89) 

The agreement of fl88D and (|89|) with the first two virial coefficients of (fl3|) in 
the large N limit is obvious for since (|88|) and (0) are exactly the same for 
any model (independently of the interaction). It is not so obvious for the second 
coefficients. First one has to note that the (constant) 0(N 2 ) term of ([18]) exactly 
cancels (|15[) . There remains the O(N) term 

ivfi + vrM-^M \ (9Q) 
V2 1 1 l + e-^H 7 v ; 

which, after Fourier transformation, indeed reproduces (|89|). 
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